NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



Si*’;. 

I*-, 

U’ 


Lockheed 
Electronics 

Company, Inc. 


ASUDSIdARyOF 
LOCKHEED CORPORATION 

1830 NASA Road 1, Houston, Texas 77058 
Tel. 713-333-5411 


InT^ : ‘ 




vof." 


lliy 


8.“®- 1 a 0-5 4/, 

NASA CR* 


Refr 644-845 
Contract NAS 9-15200 
Job Order 74-823 


TECHNICAL MEMORANDUM 

DETAILED DESCRIPTION OP THE WHEAT ACREAGE ESTIMATION 
PROCEDURE USED IN THE LARGE AREA 
CROP INVENTORY EXPERIMENT 

By 

Willa P’7. Austin 


Oiwinat photography may ha proIttsa 4 iioou 
EROS Data Center. _ 


Sioux FaW^.^ ^111 y 

comiss 

4 5 Li‘i( &1 8 8 SS 

Distribution ; 




Approved By: 


L. M. Flores, Supervisor 
Performance Analysis Section 


JSC/K. Baker 

NOAA/J. Hill 

R. M. Bizzell 

D, McCrary 

J . H . Carney 
J. L. Dragg 
R. L. Eason 

USDA/J. Murphy 

LEC/J. G. Baron 

R. B. Erb 

B. L. Carroll 

P. G. Hall 

J. E. Davis 

Ci H. Hallum 

M. L. Mannen 

P. J. Herbert 

T. c. Minter 

R. P. Ileydorn 

B. H. Moore 

R. 0. Hill 

D. G. Saile 

A. G. Houston 

P. W. Solomon 

R. D. Juday 

W. G. Straight 

G. L. Kraus 

R. E. Tokerud 

) H. L. Prior 

Technical Libri 

J. M. Sulester 

Research Data . 

D. R. Thomp.son 
M. C. Trichel 
L. C. Wade 

Job Order Pile 


February 1978 


LEC-11497 


(E80-10051) LARGE AREA CPCP IN7EHTORY N80-16399 

EXPERTHENT (LACIE) . DETAILED DESCRIPTION OE 

THE WHEAT ACREAGE ESTIHATICN PROCEDDHE HSED f*C Ar /hn ^ 

IN THE LARGE AREA CROP INVENTORY EXPERIMENT 

{Lockheed Electronics Co.) 43 p G3/43 00051 


CONTENTS 


Section Page 

1. INTRODUCTION. ..................... 1-1 

2. IMAGERY: COLOR REPRESENTATION OF DIGITAL VALUES. . . . 2-1 

2.1 COLOR METHOD OF PRESENTING PICTURE ELEMENT 

SPECTRAL SIGNATURE ........... 2-1 

2.2 FULL-FRAME IMAGERY ........ 2-2 

2.3 LACIE SAMPLE-SEGMENT IMAGERY 2-5 

3. REFERENCE MATERIALS FOR ANALYSIS 3-1 

3.1 ANALYST PACKET FOR SPECIFIC SEGMENT . ....... 3-1 

3.2 GENERAL REFERENCE MATERIAL 3-6 

4. ANALYST DECISION PROCESS. 4-1 

4.1 CHOICE OF ACQUISITIONS TO BE USED IN MACHINE 

PROCESSING . ; 4-1 

4.2 IDENTIFICATION OF SPECTRAL SIGNATURES IN SCENE . . 4-3 

4.3 DEFINITION OF TRAINING SAMPLES . . . . .... . . 4-6 

5. MACHINE CLASSIFICATION 5-1 

5.1 CLASSIFICATION 5-1 

5.2 CLASSIFICATION RESULTS ......... 5-1 

5.3 BIAS CORRECTION 5-8 

5.4 VERIFIED ACREAGE ESTIMATION SUBMITTED FOR 

AGGREGATION 5-8 

6. SUMMARY ....... ......... 6-1 

7. REFEMNCES. . 7-1 


FIGURES 


Figure Page 

1. Diagram showing analyst-satellite data interface . , . 2-3 

2. Landsat full-frame image with outline of LACIE 

segment 1528, Blaine County, Montana ......... 2-4 

3. LACIE sample-segment imagery; segment 1528, 

July 3, 1977 

{a) Product 1; channels 1, 2, and 4» . 2-6 

(b) Product 2: channels 2, 3, and 4 2-6 

(c) Product 3: channels 1, 2, and 4. 2-6 

4. LACIE sample-segment imagery: segment 1807, 

Bon Homme County, South Dakota, Jure 23, 1977 

(a) Product 1: channels 1, 2, and 4 2-8 

(b) Product 3: channels 1, 2, and 4. . . . . . . . . 2-8 

5. Flow diagram of LACIE procedure for estimating 

acreage, yield, and production 3-2 

6. Schematic diagram of detailed LACIE acreage 

estimation procedure 3-3 

7. Nominal and adjusted nominal crop calendars for 
CRD 20 (Montana) 

(a) Nominal (historical) crop calendar 3-4 

(b) Nominal crop calendar adjusted for crop 

year 1976-77 (segment 1528) 3-5 

8. Example of crop calendar adjustment from Weekly 

Meteorological Sununary 3-7 

9. LACIE sample-segment imagery for segment 1528 

at various development stages during the 1976—77 

crop year. ...... 4-2 

10. Reflectance curve for vegetation and expected spectral 
signatures for small grains 

(a) Reflectance curve of healthy green vegetation . . 4-5 

(b) Expected spectral signatures for small grains, 

segment 1528. . . . . . . ..... 4-5 


iii 


Figure 


Page 


¥ 


11. LACIE 16-channel classification products 
(segment 1528) 


(a) Plain cluster map .......... 4-8 

(b) Categorized cluster map ........ 4-8 

(c) Classification map 4-8 


12. Spectral aids, segment 1528 

(a) Scatter plots of analyst-labeled dots ..... 5-3 


(b) Trajectory plots of typical winter-wheat, 

spring-wheat, and nonwheat pixels . ....... 5-6 

(c) Green number and brightness numerical values, 

as graphed in (b) 5-7 


I 


-i} 




XV 


1 . INTRODUCTION 


During the 3 years that the Large Area Crop Inventory Experiment 
(LACIE) has been in operation, LACIE personnel have developed, 
tested, and improved the technology and procedures for estimating 
wheat acreage using data from the Land Satellites (Landsats). In 
Phase 1, techniques were evaluated for the U.S. Great Plains. 
Phase II was expanded to include estimates of acreage, yield, and 
production in the United States, Canada, and a wheat indicator 
region in the U.S.S.R. The Phase III efforts were concentrated 
on evaluating improved methods representing an increased range of 
geographic cenditions from those of Phases I and II. 

The use of satellite data to perform crop inventories utilizing 
current LACIE technology is dependent on the recognition of 
satellite-recorded radiance value patterns by an analyst. The 
subsequent machine classification and wheat acreage estimation 
accuracy are critically dependent on the analyst's correct iden- 
tification of wheat signatures. This technical memorandum docu- 
ments this basic element, the analyst-satellite interface, by 
presenting a detailed description of the LACIE process for iden- 
tifying wheat acreage. 



2. IMAGERY; COLOR REPRESENTATION OF DIGITAL VALUES 


2.1 COLOR METHOD OF PRESENTING PICTURE ELEMENT SPECTRAL 
SIGNATURE 

The Landsat-2 gathers data using an onboard multispectral scanner 
(MSS ) , which records reflected radiance values in each of four 
bands (ref. 1) ; 


MSS channel 

1 

2 

3 

4 


Wavelength ^ ym 

0.5 to 0.6 
0.6 to 0.7 
0.7 to 0.8 
0.8 to 1.1 


Spectrum 

Visible green 
Visible red 
Infrared 
Infrared 


The basic area resolution element is a picture element (pixel ) , 
which corresponds to a 0.4 5 -hectare (1.1 -acre) area on the ground. 
Thus, four values# each of which represents the average radiance 
value over 1.1 acres measured in a specific band# are recorded 
for each pixel. 


The characteristic reflectance of a pixel in the four bands is 
called a spectral signature. Spectral signatures of crops and 
other scene contents are not necessarily unique# nor are they 
constant. Signatures can change with place and time. For this 
reason# the signature identification process has not been auto- 
mated operationally and currently is done visually by analysts. 
Analysts can relate spatial and temporal information to the digi- 
tal radiance values for signature identification. To facilitate 
this analyst identification application# digital values are 
represented as colors and presented in imagery form. Three 
colors are used (red# blue, and green) , so information from three 
channels or bands can be presented at one time. Color represen- 
tation of digital values is not perfect# but it is very effective 
as a vehicle for providing information. 



A 


MS3 data aro telemetered to the National Aeronautics and Space 
Administration (NASA ) , Goddard Space Flight Center (GSFC) . The 
GSFC forwards the satellite data on film to the U.S. Department 
of Agriculture (USDA) and on computer-compatible tape (CCT) to 
the Lyndon B. Johnson Space Center (JSC). Figure 1 summarizes 
the analyst-satellite data interface. 

2 • 2 full-frame imagery 

A Landsat full-frame image corresponds to a 185- by 185-kilometer 
(100- by 100-nautical-mile) area on the ground. Full-frame 
imagery is created at the USDA Aerial Photogrammetry Field Office 
at Salt Lake City, Utah, from 70-millimeter archival film supplied 
by GSFC. Figure 2 is a full-frame image of data acquired July 3, 
1977, over Montana. Channels 1, 2, and 4 are used in the com- 
posed false-color image; and channel color assignments are blue, 
green, and red, respectively. 

The satellite telemeters the radiance values to GSFC in the 
range 0 to 63 for each channel. At GSFC, channels 1, 2, and 3 
are decompressed logarithmically into the range 0 to 127; 
channel 4 does not require decompression. Values are mapped 
directly onto a film positive calibrated to a 16-level-density 
gray scale: 0 into black, and 127 (or 63) into white. The full- 

frame color composite is created from these black-and-white images 
of the individual channels. Full-frame imagery for every Landsat 
acquisition containing a LACIE segment is provided to JSC for 
analyst use at least once per biowindow; to fulfill this require- 
ment, cloud cover must be less than 20 percent on the imagery. 

This imagery is used to monitor episodic events and to give the 
analyst a large view of the geographical area which includes a 
LACIE segment. Full-frame imagery usually clarifies the distinc- 
tion between agricultural and nonagricultural areas. 



YIELD AND PRODUCTION ESTIMATES. 


Figure 1.— Diagram showing analyst-satellite interface 



NASA-S-77-12394 



.andsat tull-frame iir.age with outline of LACIE 
egment 1528, Blaine County, Montana. 













2 . 3 LACIE SAMPLE-SEGMENT IMAGERY 


A 9- by 11-kilometer (5- by 6-nautical-raile) LACIE sample segment 
(segment 1528 in Blaine County^ Montana) is outlined within the 
full frame illustrated in figure 2. 

At GSPC^ an 18.5- by 20-kilometer (10- by 11-nautical-mile) search 
area, based on the geographical coordinates of the segment, is 
extracted from the full-frame digital data and used as a basis 
for an image edge detection routine to locate exactly and extract 
the 9- by 11-kilometer LACIE sample segment (ref. 2). In the 
first process, bias and scale factors are computed from the mean 
value and standard deviation of each channel over the entire 
search area. Then, the bias and scale factors are placed in the 
headers of the CCT's which contain the four digital values for 
each of the 22 932 pixels in the sample segment; the digital 
values are extracted from the second process. These CCT's are 
sent to JSC. Sample segment imagery, which is a color representa- 
tion of these digital values, is created at JSC on the production 
film converter (PPC) using the bias and scale factors. 

Figure 3(a) is segment 1528 outlined on figure 2. This film 
product was created from digital values in channels 1, 2, and 4 
(two channels in the visible range and one in the infrared range); 
color assignments, as in figure 2, are blue, green, and red, 
respectively. This type of image is termed LACIE product 1. 

A film product created with channels 2, 3, and 4 (two channels 
in the infrared range and one in tne visible range) is made so 
that information in channel 3 is available for the analyst . 

Color assignments in this product are red (channel 2), blue 
(channel 3) , and green (channel 4); polarities are reversed for 
channels 3 and 4. This product, LACIE product 2, is shown in 
figure 3(b). 
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(a) Product 1: channels 1, 

2, and 4. 


(b) Product 2: channels 2 , 

2 , and 4. 




(c) Product 3: channels 1, 

2 , and 4 . 
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Figure 3.— LACIF sample-segment imagery: segment 1528, 

July 3, 1977. 
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The analyst uses imagery for spatial information and for spectral 
value information. The film products shown in figures 3(a) and 
(b) are generated to emphasize contrast. These products are 
excellent for field delineation and enhanced spatial features. 
However, depending upon the data in the scene, contrast is some- 
times achieved at the expense of consistent color depiction of 
spectral values. Figure 3(c), LACIE product 3, is imagery devel- 
oped specificaily for more consistent color display of spectral 
signatures; this is essential for accurate identification. Fig- 
ures 4(a) and (b) , imagery of segment 1807, provide an additional 
comparison of LACIE products 1 and 3 . The channels used and the 
color assignments for product 3 are the same as for product 1, but 
the method of depicting digital values as colors is different. 

For product 1 [figs. 3(a) and 4(a)] and product 2 tfig, 3(b)], 
each channel is mapped individually onto the designated color- 
density scale. The bias and scale factors computed at GSFG are 
used to expand the data range for each channel over the entire 
color range. Interchannel relation is not necessarily preserved. 
For product 3 [figs. 3(c) and 4(b)), an overall mean value for 
the three channels is computed (channel 4 values are doubled 
before computing this) . The overall mean is set at the midpoint 
of the color-density scale. No bias is used, but there is a scale 
factor. The entire range of color is not necessarily used for each 
channel, and interchannel relationship is preserved (refs. 3, 4). 

Figures 3(a) and (c) are very similar; however, figures 4(a) and 
(b) are markedly different. In this case, maximizing contrast 
has altered the relationship between the channels and changed 
the spectral signature representation of many of the pixels. 

Figure 4(e) would be very misleading if used for signature 
identification. 

The use of imagery as a means of implementing the analyst- 
satellite data interface is fundamental in LACIE. There are 
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(a) Product 1: channels 1, 

2 , and 4 . 



(b) 
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Product jt channels 1, 
2 , and 4 . 
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Fiqurc 4.— LACTF sample-seqmont imaqory: seqment 

Don Homme County, South Dakota, June 2 ^, 197o 
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91 728 digital values for each LACXE sample segment (one value 
in each of four bands for each of 22 932 pixel locations) . If 
presented in numerical form, this would be overwhelming. For 
each Landsat acquisition# the analyst has available the 9- by 
11-kilometer segment imagery; that is# product 1 (channels 1, 2, 
and 4)# product 2 (channels 2# 3# and 4)# and product 3 (chan- 
nels 1# 2, and 4). Pull-frame imagery coverage is usually avail- 
able for the previous year's acquisitions and occasionally for 
the current acquisition. Imagery provides a practical means of 
identifying spectral signatures and separating spatial features. 


3. REFERENCE MATERIALS FOR ANALYSIS 


The LACIE approach for estimating acreage, yield, and production 
is shown in figure 5. The average proportion of wheat for sample 
segments within the stratum is used to obtain an estimate of the 
stratum proportion of wheat. This is multiplied by the stratum 
agricultural area to obtain an estimate of the stratum wheat 
acreage. In LACIE, strata are defined by percentages of agricul- 
tural density using Landsat data, soil information, and climatic 
conditions. A simple random sample within each stratum is used, 
and a predetermined number of sample segments of specified loca- 
tion are allocated per stratum, 

A more detailed view of the exact procedure used to produce an 
acreage estimate is diagrammed in figure 6. With this diagram 
as a basis, the following subsections present a review of mate- 
rials available to aid tiie analyst in segment analysis. 

3.1 ANALYST PACKET FOR SPECIFIC SEGMENT 

LACIE personnel at JSC have assembled an analyst packet for every 
sample segment. The full-frame imagery is used to establish the 
relationship of the segment to the area specified in the reference 
material contained in the packet. 

The packet contains the following data. 

a. Maps: 1:24 000 and 1:250 000 scale. 

b. Ancillary data: 

• Cropping practices for the Crop Reporting District (CRD) . 

• Soils data for the CRD. 

• Nominal crop calendar for the CRD based on 10-year averages 
[fig. 7(a)]. Information on average length of crop develop- 
ment stages, the nominal dates when these stages occur, and 



Figure 5.— Flow diagram of LACIE procedure for estimating acreage, yield, and production. 
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PLT • planting 
EMG ■ emargance 
JNT • lointing 
MED ■ heading 
SD • toft dough 

TRN • turning 
R ■ ripe 

MAT ■ maturity 
HRV ■ harvett 

E under stage ■ rough estimate of date 





(a) 

Nominal 

(historical) crop calendar. 



Figure 7.— 

Nominal 
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and adjusted nominal crop calendars 
CRD 20 (Montana) . 
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CROP CALENDARS PLOTTED 01/06/76 
PERCENT OF AREA IN DEVELOPMENT STAGE BY SPECIFIED DATE FOR 
MONTANA RECENT AVERAGE AVERAGE CROP CALENDARS CRD 70 



R - np« 5/28 ad|utt«d to BMTS 3.5 (W) 


MAT • maturity 
HRV ■ harvtit 

E undtr itaga ■ rough ettimatt of data 


(b) Nominal crop calendar adjusted for crop 
year 1976-77 (segment 1528) . 


Figure 7.— Concluded. 
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the relative growth stages o£ other crops in the area is 
presented. Figure 7(b) shows adjustments to the nominal 
crop calendar. By correlating the year's specific growth 
stages found on the crop calendar adjustment (fig* 8) from 
the Weekly Meteorological Summary (ref. 5 ) , the analyst 
adjusts the nominal crop calendar to be more specific for 
the current year, as shown in figure 7(b). 

• Historical crop percentages for the political subdivision 
for the previous 4 or 5 years. 

c. Imagery for the segment: 

• Product 1 for all acquisitions during the past crop year, 
if available. 

• Products 1, 2, and 3 for each acquisition in the current 
crop year. 

d. Machine classification data: 

• Previous wheat acreage estimates on the segment. 

• Classification products for the most recent estimate, 

3.2 GENERAL REFERENCE MATERIAL 

Reference materials not included in the analyst packet are 

available: 

a. Weekly Meteorological Summary (ref. 5) — These summaries pro- 
vide current data on temperature and precipitation, a crop 
calendar adjustment which reflects the current Robertson 
Biometeorological Time Scale (BMTS) growth stage of wheat, 
and a narrative summary for crop and weather assessment on a 
statewide basis. Episodic events, such as recent rainfall, 
alter crop signatures; soil reflectivity contributes an 
indeterminate component to the average reflectance value for 
an acre recorded by Landsat, and the reflectivity of wet soil 
differs from that of dry soil. Long-term events, such as 



Tho kiootniiing of each growth itago Is when approximately 50 percent of the crop has reoched one of the following numbered stages, 

Growth stage Des cription Growth stage 

1.0 Planting * 4.0 Heading 

2.0 emergence ^ 5,0 Soft dough . ' ^ ? Vfi' TH R 

Dormancy (winter wheat) 6,0 Ripe iiui .ij. ■* i ^ > . 

Sprlnggrowth (winter wheat) 7.0 Harvest PAGIj X > POOR 

3.0 Jointing (greenlng'Up period) wavav* 


Figure 8.- Example of crop calendar adjustment from Weekly 
Meteorological Summary. The numbers refer to the 
Robert son BMTS growth stages for wheat. 


drought, affect the expected Robertson BMTS for wheat, as 
well as expected spectral signatures. 

b. Landsat full-frame Imagery. 

c. LACIE Phase III Classification and Mensuration Subsystem 
(CAMS) Detailed Analysis Procedures (ref. 6) — This guide 
is used by all analysts to ensure a consistent decision 
process . 

d. Analyst Interpretation Keys (ref. 7) - A definitive reference 
for crop signatures, vthio document gives a representative set 
of segments over a stratum, along with crop calendars, full- 
frame imagery, and labels for selected fields. 


If 


4. ANALYST DECISION PROCESS 


The analyst correlates the available material. Two decisions 
must then be made prior to machine classification of the segment; 
(1) selection of the acquisitions to be used for machine process- 
ing and (2) specification of the training data to be used as the 
basis for machine classification. Reference materials are usually 
consulted at many stages in the decision process. Decisions may 
be updated as frequently as every 18 days, corresponding to each 
successive jjandsat overflight. 

4.1 CHOICE OF ACQUISITIONS TO BE USED IN MACHINE PROCESSING 

The analyst may choose to classify a segment using a single acqui- 
sition (four channels of data) . however, multitemporal classifi- 
cation (classification using data frofa more than One acquisition) 
is preferred because it utilizes the temporal separability of the 
growth cycle of wheat from the growth cycles of confusion vege- 
tation. A maximum of 16 channels of data, or 4 acquisitions, 
can be used for machine classification. 

Acquisitions to be machine processed are chosen according to the 
following criteria. 

a. Imagery quality. Cloudy or hazy imagery is omitted if pos- 
sible. Acquisitions used for machine classification must be 
mutually registered. 

b. Wheat growth stages. 

c. Wheat separability from confusion crops should be maximized 
by the acquisitions selected. 

For example, figure 9 shows a sequence of acquisitions for seg- 
ment 1528, four of wnicn are selected for classification; fig- 
ure 2 is the full-frame imagery for segment 1528, acquisition 
date July 3, 1977; and tne crop calendars for the segment are 
illustrated in figures 7(a) and (b) . This segment is a 
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November 11, 1976 December 17, 1976 



riyure 9.— Sample-segment imagery for segment 1528 at various 
development stages during the 1976-77 crop year. 
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mixed-wheat segment, so both winter- and spring-wheat growth 
patterns must be considered in choosing acquisitions. 

According to the crop calendar, winter grain was the only crop 
in a growth stage on November 11, Stream growth is also evident 
on the imagery, but multitemporal processing (used wherever 
possible) should eliminate this confusion factor. According to 
the crop calendar, winter wheat was in the dormancy stage on 
December 17. Based on the acquisition-selection criteria, either 
the November 11 or the December 17 acquisition would be a logical 
choice for processing. The April 22 imagery was acquired at both 
the planting stage of spring wheat and a growth stage of winter 
wheat. This also would be a good processing choice. Notice that 
if multitemporal processing were done on two of these dates alone 
(for example, November 11 and April 22) winter grains would be 
confused with stream growth. The streams would have to be desig- 
nated other (DO) , which is an area outlined to be designated as 
nonwheat, regardless of the statistical classification results 
within the area. 

The May 28 and July 3 acquisitions are both of good image quality 
and add separability from natural vegetation (the full-frame 
imagery shows clear separation) and information on the growth 
stage of wheat. One of these dates should be used. 

On August 8, the wheat was at a critical stage in development; 
thus, this acquisition would be a logical choice for inclusion 
in multitemporal machine processing. At that time, small grains 
(winter and spring) were ready for harvest or had been harvested, 
and confusion vegetation could be separated from small grains. 

4.2 IDENTIFICATION OF SPECTRAL SIGNATURES IN SCENE 

Using all available information, the analyst identifies the spec- 
tral signatures in the segment. This is done by correlating the 


temporal growth stages o£ small grains and the expected sequence 
of changes in e;;iectral signatures to the sequence of spectral 
signatures evident on the imagery. Plant phenology is indicated 
by a temporal change in the infrared reflectance relative to the 
reflectance in the visible range, A growing plant has relatively 
high reflectance in the infrared channels because of scattering 
from the internal structure of the leaf; green plants have rela- 
tively low reflectance in the visible range (green, red) because 
of absorption by chlorophyll (ref, 8), Figure 10(a) shows a 
typical curve for green growing vegetation plotted in wavelength 
versus percentage of reflectance, Since channel 4 (0,8 to 
1.1 micrometers) is assigned red in LACIE products 1 and 3 
imagery, a red spectral signature on this imagery will reflect 
healthy green vegetation. Chlorophyll is sensitive to changes 
in plant condition, such as stress and senescence, lower chloro- 
phyll content in vegetation results in less absorption in the 
visible range, thus giving more relative reflectance resulting in 
a less red signature. Figure 10(b) shows the expected sequence 
of small-grain signatures for the acquisition dates of figure 9, 

For example, follow the field labeled W on the sequence of imagery 
in figure 9. This is clearly a winter-grain field. The field 
labeled S follows a spring -grain temporal growth pattern. The 
fields labeled N do not follow the temporal grfenfth pattern of 
wheat . 

Any events which might affect the evident signatures, such as 
drought, recent rain, and the irrigated mode of crop practices, 
must be allowed for in the expected signature. 
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12 3 4 WAVEIENOTH (MICRONS) 

(a) Reflectance curve of healthy green vegetation. 


DATE 

WINTER GRAIN 

SPRING YCHEAT 

NOV. t1 10 
MAR. 30 

PLANTING/EMERGENT STAGE - 
NO RED TO RED 

PREPLANTING - NO RED 

APR. 22 

VIGOROUS GROWTH STAGE - RED 

PLANTING -NO RED 

MAY 28 

VIGOROUS GROWTH STAGE 
perhaps SOME CHANGE 
PRcPARATORY TO TURNING ~ RED 
OR BRICK RED 

EMERGENT- RED OR PINK, 
DEPENDING ON PLANT 
CANOPY 

JULY 3 

RIPE STAGE - ORANGE. YELLOW, 
BROWN 

VIGOROUS GROWTH STAGE - 
RED 

AUG 8 

HARVEST STAGE - YELLOW, WHITE, 
TAN 

HARVEST OR READY FOR 
HARVEST - YELLOW, BROWN, 
OLIVE GREEN, WHITE 


(b) EJcpected spectral signatures for small grains, segment 1528. 


Figure 10.— Reflectance curve for vegetation and expected 
spectral signatures for small grains. 




4 . 3 DEFINITION OF TRAINING SAMPLES 


When signatures have been identified and acquisitions have been 
chosen for machine processing (acquisitions of good quality, 
which are registered and representative of the wheat growth 
cycle, as well as capable of gaining maximum separability of 
wheat from confusion factors in the scene) , training samples 
must be defined tor machine classification. 

A subset of the 22 932 pixels in the LACIE sample segment is 
identified and used to train the classifier. In crop year 
1975-76, the analysts chose training fields from the spectral 
signatures in the scene? analysts used training fields to define 
the stratification of Landsat spectral space. Several fields, 
preferably in different geographical areas of the sample segment, 
were selected to approximate a spectral signature. Fields 
approximating the same spectral signature were grouped by the 
analyst into subclasses. The subclass mean value and standard 
deviation in each channel were calculated, and all pixels in the 
segment were classified using these subclass statistics generated 
by the training fields. 

For crop year 1976-77, the analysts identified individual pixels; 
spectral stratification was done using a machine clustering 
process. The image display of a LACIE sample segment has 
196 pixels per line and 117 lines; a grid mark appears at the 
intersection of every 10 pixels and 10 lines. The pixel at the 
lower right corner of each block of 100 pixels so delineated has 
been designated a dot. Two preselected random samples of dots 
were specified to be labeled by the analyst: Type 1 dots were 

used to initiate the clustering process and to label the clusters 
generated, and type 2 dots were used for bias correction of the 
final wheat acreage estimate (section 5.3). The analyst labels 
a minimum of 30 dots from the random selection specified for 
type 1 dots. These dots (as in training field selection) were 




selected to be good samples of spectral signatures; no border or 
edge pixel^/ were used. Twenty of these labeled dots were selected 
randomly to initiate clustering. 


Clustering is the process of grouping pixels according to some 
distance measure. The pixel vector (X = Xj^,X 2 , • • * >x^; n = number 
of channels) of each of the 22 932 pixels in the segment was com- 
pared with the pixel vector of each of the 20 starting dots 
(Y = yj^/y 2 » • * * • Each pixel is assigned to the closest start- 

ing dot. 


Distance = 




where m is the number of channels used in clustering. 


After all pixels have been assigned, the mean and standard devia- 
tion of each cluster are computed for each channel. Clusters can 
"split and combine" using predefined parameters, to a maximum of 
60 clusters. A maximum of 10 passes can be made through the data, 
assigning pixels to clusters, recomputing cluster statistics, and 
then splitting clusters (if channel standard deviations are 
large) , merging clusters (if intercluster distances are small) , 
or deleting clusters by eliminating cluster mean vectors (if the 
cluster populations are too small). After the first pass, pixel 
vectors are "distanced" from cluster mean vectors (refs. 9, 10). 
Figure 11(a) is a color-coded cluster map of segment 1528. 

Figure 11(b) is a categorized (conditional) cluster map. The 
mean vector of each cluster has been distance-compared to the 
dot vector of each of the 30 analyst-labeled type 1 dots and 
labeled wheat or nonwheat according to the label of the closest 
dot. If the distance was within a specified threshold of a dot 
labeled S, W, or N (spring wheat, winter wheat, or nonwheat), the 
cluster was color-coded green, cyan, or yellow, respectively. 





If it was not within the threshold, category assignment of the 
cluster was according to the closest dot, but color coding was 
unique. The cluster mean value and standard deviation in each 
channel, with cluster labels assigned as above, are used to 
classify the segment. 


5. MACHINE CLASSIFICATION 


5.1 CLASSIFICATION 

Pixel-by-pixel classification of all pixels in the segment is 
performed using the LACIE System (ref. 9), Either subclass or 
cluster statistics are used to build the LACIE classifier for a 
segment. 

A classifier is a system of subclass density functions, one 
function for each cluster or subclass. The probability of each 
pixel in the segment belonging to each subclass is computed; 
that is, the pixel vector is substituted in each subclass density 
function and the value of the function is calculated. Values are 
summed to the category level. Pixel assignment is made to the 
most likely category and then to the most likely cluster or sub- 
class within that category. Pixel assignment is not directly to 
the most probable subclass, since the CAMS procedures specify a 
summation-of-likelihoods classifier. Classification is on a 
pixel-by-pixel basis; that is, the procedure above is repeated 
22 932 times in order to classify a segment; all pixels are 
classified. 

After classification, thresholding is applied to remove from 
final classification results pixels with low probability of 
belonging to the assigned subclass (ref, 11) . 

5 . 2 CLii.SSIFICATIQN RESULTS 

Classification results are presented in the forms of imagery, as 
shown in figure 11(c), and numerical results. The imagery repre- 
sentation permits the analyst to do a spatial verification (clean 
field patterns; wheatfields classified as wheat; and natural 
vegetation, streams, and other crops classified nonwheat) . If 
areas are misclassified, the cluster maps [figs. 11(a) and (b) ] 
permit monitoring of the machine clustering routine and hence of 


the segment classifier, if classification was based on cluster 
statistics. On the figure 11(c) classification map, light gray 
is used for pixels classified nonwheat, medium gray for pixels 
classified spring wheat, and black for pixels classified winter 
wheat. Pixels that were thresholded after classification are 
white on this classification map. Numerical results, wheat 
acreage estiraation for the segment and the classified category 
of each of the dots or training fields, according to the proce- 
dure used to select training samples, are presented. This permits 
further verification of results in accordance with well-defined 
procedures (ref. 6). 

Spectral aids, based on a transformation of the data (refs. 12, 
13), are provided with classification results. These enable the 
analyst to check the consistency of dot labeling. Spectral data 
are transformed into two variables: green number and brightness. 

This information is presented to the analyst in the following 
forms. 

a. Scatter plots - graphic representation of data distribution 
and the distribution of analyst dot labels for each acquisi- 
tion used in the classification. Figure 12(a) shows scatter 
plots of analyst-labeled dots for segment 1528, acquisitions 
of April 22, 1977; May 28, 1977; July 3, 1977; and August 8, 
1977. 

b. Trajectory plots - time traces of an individual pixel. One 
trajectory plot is generated for each of the 209 dots when- 
ever multitemporal processing is used. Figure 12 (b) shows 
sample trajectory plots for a winter-wheat, a spring -wheat, 
and a nonwheat dot. 

c. Numerical Values for the dots used as representative trajec- 
tory plots in figure 12(b) are presented in figure 12(c). 
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(a) Scatter plots of analyst-labeled dots. 
Figure 12.- Spectral aids, segment 1528. 
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Figure 12 Continued, 
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Figure 12.— Continued. 
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Acquisition 1; 
Acquisition 2: 
Acquisition 3: 
Acquisition 4: 


Apr. 22, 1977; BMTS 2.5, winter wheat, and 2.0, spring wheat 
May 28, 1977; BMTS 3.5, winter wheats and 3.1, spring wheat 
July 3, 1977; BMTS 5.0, winter wheat, and 4.6, spring wheat 
Aug. 8, 1977; BMTS 7.0, winter wheat, and 6.0, spring wheat 


(b) Trajectory plots of typical winter-wheat, spring-wheat, 
and nonwheat pixels. 


Figure 12.— Continued. 


DOT 

NUHDER 

GRID 

INTERSECTIONS 

analyst 

lASEL 

CLASSiriED 

ACQUISITION 1 
G R 

ACQUISITION 2 ACQUISITION 3 ACQUISITION 4 
G B 6 D G B 

91 

50 150 

W 

w 

5 

65 

13 

74 5 

80 5 93 

119 

70 50 

S 

S 

Z 

61 

11 

86 28 

72 6 69 

1S3 

90 10 

N 

N 

2 

66 

2B 

65 18 

65 22 58 


(c) Green number and brightness numerical values, as 
graphed in (b) . 

Figure 12.— Concluded. 



5,3 BIAS CORtUSCTION 


A classifier bias is introduced in the classification machine- 
decision process. In crop yeaL 1976-77, a statistical correction 
factor was applied to the wheat acreage estimation generated by 
machine classification (ref, 14). To correct for the classifier 
bias, an independent, preselected, random sample of type 2 dots 
was labeled by the analyst before machine classification (sec- 
tion 4.3). Each dot specified for this random sample had to be 
labeled, including border and edge pixels; 60 random dots were 
labeled for bias correction if one category of wheat was used; 

99 dots were labeled if both spring and winter wheat were present. 
Acreage estimates as classified and as bias corrected are pre- 
sented to the analyst for verification. 

5.4 VERIFIED ACREAGE ESTIMATION SUBMITTED FOR AGGREGATION 

An independent check of all dot labels to ensure consistency in 
labeling was made by a quality assurance team of analysts which 
did not include the analyst to whom the segment was assigned. If 
verification of any area by the analyst or by the quality assur- 
ance team indicated that a rework would improve accuracy of 
results, this was done. 

The verified acreage estimation for the wheat in the sample seg- 
ment (in crop year 1976-77, this was a verified bias-corrected 
acreage estimate) is submitted for aggregation and subsequent use 
in production estimates (fig. 5) . 


r 




6 . SUMMARY 

The use of satellite data fot crop inventories, as in LACIE, cur- 
rently demands intense interaction of the human analyst with the 
digital data from the satellite. This is implemented by usihg 
imagery forms of the digital values. Analysts must be skilled 
in image interpretation, as well as in the machine processing 
techniques applied to the digital values. 

Data acquired by the Landsat MSS are telemetered to GSFC, where 
LACIE sample-segment digital information is extracted. An ana- 
lyst at JSC correlates a variety of ancillary data and applies 
this to an imagery representation of the digital data for identi- 
fication of sample-segment scene content. The analyst labels a 
subset of the 22 932 pixels in the sample segment. Based on 
this subset of labeled pixels, machine classification of all 
22 932 pixels is performed using the digital values. The esti- 
mate of wheat acreage proportion for the sample segment is com- 
puted. This sample-segment proportion estimate is used to 
estimate wheat acreage proportions over larger geographical 
areas and to produce estimations of production. 

Experience gained from the 3 years of development in LACIE can 
serve as a foundation for further progress in the use of satel- 
lite data. A detailed examination of the LACIE analyst procedure 
for identifying wheat and estimating wheat acreage can provide a 
basis for redefining the analyst-satellite interface, perhaps 
toward more automated, less intense interaction. 
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